optimal laser pulse design for transferring the coherent nuclear wave packet of 
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Within the Franck-Condon approximation, the single ionization of H2 leaves H.^ in a coherent 
superposition of 19 nuclear vibrational states. We numerically design an optimal laser pulse train to 
transfer such a coherent nuclear wave packet to the ground vibrational state of H J . The simulation 
results show that the population of the ground state after the transfer is more than 91%. Frequency 
analysis of the designed optimal pulse reveals that the transfer principle is mainly an anti-Stokes 
transition, i.e. the in lsa g with excited nuclear vibrational states is first pumped to 2po g state 
by the pulse at an appropriate time, and then dumped back to lsa g with lower excited or ground 
vibrational states. 

PACS numbers: 33.80.Rv, 42.50.Hz, 02.30.Yy 



I. INTRODUCTION 



Controlling coherent quantum states has been a long- 
standing goal since the invention of laser pulses. With the 
rapid advent of technology in recent years [l[ , researchers 
can now fine-tune the laser parameters to control the ul- 
trafast processes inside atoms and molecules For 
example, by varying the relative phase of a two-color (w- 
3w) laser field with w the fundamental angular frequency, 
the target molecule may be constructively or destruc- 
tively excited by simultaneously absorbing the u and 3w 
photons Q . Another example is to change the time delay 
between two laser pulses so that the molecule can be first 
pumped to an intermediate state, then evolves and accu- 
mulates the time-dependent phases, and later be dumped 
to a different final state, thereby changing the production 
of a chemical reaction |6j|. Furthermore, specific tailoring 
of the laser field may dictate a complex chemical reac- 
tion to follow one particular channel and stay away from 
all the others, achieving a selective terminal state Q. 
Most recently, thanks to the phase-stabilized few-cycle 
laser pulse §, the emission of an ionized electron J{|, 
or the charge-direct transfer between nuclei [13, HH also 
become possible. The attosecond pulse may be used to 
selectively excite or ionize the target at unprecedentedly 
precise timing during the chemical reactions (i~2l-[l4| , and 
thus helps understanding the time-resolved fundamental 
physics. 

As the simplest neutral molecule, H2 (or D2) is often 
chosen as a prototype system to be controlled and ana- 
lyzed. In the past few decades, the basic processes for 
H2 exposed in strong laser fields have been extensively 
studied. As shown in Fig. [TJ after one electron absorbs 
enough photon energy and escapes from the nuclei, it is 
left with a molecular ion Ht in lsa„ state UM. From 




FIG. 1: (Color online) The schematic picture for the interac- 
tion between laser pulses and H2. Four curves from bottom up 



are the potential curves for H2, Hj 



in 2pa u , and 



coulomb explosion. F-C remarks the Franck-Condon transi- 
tion area. The single ionization of H2 produces the in 
lsa g , followed by the dissociative ionization by the probe 
pulse. The kinetic energy release (KER) reflects the infor- 
mation when the coulomb explosion happens. 



the Franck-Condon approximation [16j . we can assume 
that initially the nuclear wave packet (NWP) of is 
the same as the ground state of H2 , and then it evolves 
along the lsa g potential curve, as shown in Fig.rjjc) [TtI ] . 
If a time-delayed probe pulse is subsequently introduced, 
Hj" may dissociate through the laser-induced coupling 
between lsa g and 2pa u [ill. Il8l - [2l| . The mixture of the 
dissociative channels, i.e. the paired and unpaired states, 
will induce an asymmetric electron localization [l2|, [22| — 
|3~0 | . Alternatively, may also be ionized by the probe 
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pulse, which leads to the Coulomb explosion [151 125 
The intcrnuclcar distance when the ionization of tak- 



2 



ing place can be reflected by the kinetic energy release 
(KER) of the Coulomb-explosion fragments |33l - [3a |. In 
addition to these non-electron correlation processes, the 
first ionized electron may come back and rescatter with 
H 2 , accompanying with the excitation of (37l l38ll and 
high harmonic generation [39l ] , or auto-ionization (40j . If 
the laser environments are appropriate, the single ioniza- 
tion of H2 may leave the H J in higher excited electronic 
states, e.g. 2pa u [l4| or 2pn u |4lj . 



In these processes, the complexity of the NWP 
makes the whole process even more complicated. After 
the single ionization, the NWP of is a superposition 
of 19 vibrational states with negligible auto-dissociative 
states. Each vibrational state has a different spatial dis- 
tribution. Once the molecule is dissociated or ionized, 
each vibrational state also contributes a different KER 
to the molecular fragments J42|. The coherent super- 
position of the vibrational states may partly smear the 
asymmetric electron localization |31| , or lead to the time- 
dependent dissociation [43| . Since the stationary nuclear 
state can significantly simplify the physical picture, it is 
often desired to transfer the coherent NWP to one sta- 
tionary state, especially to the ground vibrational state 
ofH+ 



For the hetero- nuclear molecular ion HD + , the perma- 
nent dipole induced by the asymmetric nuclear mass may 
transfer the Franck-Condon NWP to the ground vibra- 
tional state [44j|. However, for H^, the external laser field 
has to be applied. Niederhausen and Thumm suggested 
to use the multi-pulse protocol to control the Franck- 
Condon coherent NWP and found that in the final coher- 
ent population, the largest proportion for a certain vibra- 
tional state can exceed 60% |45||. Niikura et al. studied 
to exert a laser-induced dipole force at an appropriate 
time to achieve up to 50% population for the ground 
vibrational state [46|. Picon et al. proposed to use a 
chirped fcw-hundrcd-fcmtosecond pulse or pulse train to 
transfer the first and second excited vibrational states 
to the ground vibrational state with the proportion up 
to 90% (47) ■ Bryan et al. used the pump-modify- probe 
strategy to manipulate the vibrational states, where the 
time-delayed second pulse may modify the relative pop 
ulations of different states 14 



In this paper we use optimal control theory to numer- 
ically design a laser pulse train to tailor the coherent 
vibrational states. We formulate it as a minimax prob- 
lem with bounded constraints, and then apply sequen- 
tial linear programming algorithm [49| to solve it. The 
gradient of the performance metric with respect to the 
laser pulses can be derived in an analytic manner, which 
facilitates the numerical computation. For the initial 
Franck-Condon NWP, the transfer to the ground vibra- 
tional state is achieved with the population more than 
91%. 
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FIG. 2: (Color online) (a) Franck-Condon coefficients of the 
NWP of H^. (b) Wave functions for the 19 vibrational states 
(in linear scale), (c) The free propagation of the Franck- 
Condon NWP of in the ls<r 3 potential curve (in loga- 
rithmic scale). 



II. NUMERICAL MODEL 
A. Two-state equation 

Consider the case that the single ionization of H2 pro- 
duces a free electron and a molecular ion H+ in 1 
where the NWP of H is described by the Franck-Condon 
approximation. If the time-delayed probe pulse is intro- 
duced to cause the dissociation, the dynamics is mainly 
governed by a two-state equation (atomic units are used 
unless otherwise stated) 







~1pg(R,t) 









T R + V g (R) d gu {R)E{t) 
d gu {R)E{t) T R + V U (R) 



i> g (R,t) 



(1) 

where ijj g (R, t) , ip u (R, t) are the NWP corresponding to 
the electron in lsa g and 2p<r u states, and V g (R), V U (R) 
are the potential curves for lsa g and 2pa u states, respec- 
tively. The dipole coupling between these two states is 
represented by d gu , and Tr = — j|y is the second or- 
der differential operator, where M = 918 is the reduced 
mass of two nuclei. The molecular rotation is neglected 
since we limit the pulse duration within a few tens of 
femtoseconds. The initial NWP is given by 



V s (i?,o) = v£ r (#), MR0) = Q, 



(2) 



where ipQ T (R) is the ground state of H2 . We use the Split- 
Operator method to solve Eq. ([T]) [Ho| . The R spans from 
to 40, and the spatial step AR = 0.04. The time step 
is set as At = 1 . Mask functions are used to suppress the 
unphysical reflection by the boundary of the simulation 
box. 
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Our objective is to design a laser pulse E such that at 
the terminal time Tf , ip g can be transferred to the ground 
state of Ho", i.e. 



(3) 



where v is the index of the vibrational state. The ini- 
tial NWP ip g (R, 0) is mainly a superposition of 19 vibra- 
tional states. By projecting it to the vibrational eigen- 
states of Hj, we obtain the Franck-Condon coefficients, 
as shown in Fig.[2]Ja). The three vibrational states v = 1, 
2, and 3 amount to around 50% of the total population. 
Fig. [U[b) shows the wave function for all these 19 vibra- 
tional states, and Fig. [^c) plots \ip g (R, t)\ 2 . The black 
curve in Fig. [^Jc) is the expected time-dependent inter- 
nuclear distance (R(t)). Clearly, the NWP goes through 
a collapse and revival procedure, and the revival time is 
about 300 fs [HG3. 



B. Numerical optimization algorithm 

We formulate the design of a laser field E to realize 
the NWP transfer of tp g as a minimax problem, and then 
apply a sequential linear programming algorithm to solve 
it. To avoid ionization of H^, we restrict the amplitude 
of E within 0.1, and the pulse duration less than 32 fs. 

The wave function transfer is formulated as a con- 
strained minimax problem on the laser electric field E: 



subject to 



min max J n , 

E n.e{0,l,-" ,iV-l} 



\E(t k )\ <0.1, k = 0,--- ,K-l, 



where 



J n 



^°JR n )~e- ia iP g {R n ,T f )\ 



(4) 



(5) 



(6) 



Here n (or k) is the index for the spatial (or temporal) 
step, N (or K) is the total points in the spatial (or tem- 
poral) axis, and a is a global phase to be determined 
soon. The function J n quantifies the difference between 
the desired and actually achieved wave functions at the 
spatial grid R n . If the maximum error of J n is mini- 
mized over the whole spatial range, one can expect that 
the desired wave function is achieved. 

The NWP transfer fidelity can be measured by 



F = Re{e ia {i> g {R,T f W=\R))} 



(7) 



The global phase a in Eqs. (J6j) and (O can be obtained 
by maximizing the fidelity: 



arg 



{(TP g {R,T f W=\R))}, 



(8) 



where arg denotes the argument of a complex number. 

The minimax problem has been extensively studied in 
the optimization and control community [5ll453l | . To find 



the optimal laser field, we start from an initial guess and 
then gradually approach the optimal solution by itera- 
tion. Suppose that at the j-th iteration, the current laser 
pulse is E 3 . We need to determine a small increment 
AE 3 such that at the (j + l)-th step, the new laser pulse 
E 3+1 = E 3 + AE 3 is a better solution to minimize the 
transfer error J n . By first order approximation, we have 



Jn(E 



J n (E 3 )+V T E] J n (E 3 )AE 3 . 



(9) 



The analytic derivation of the gradient \7'^ j J n (E 3 ) is 
given in the appendix. 

We then apply a sequential linear programming algo- 
rithm as follows: 

1. Choose a small constant as the initial guess of the 
electric field; 

2. At the j-th step, compute J n (E 3 ) and V^j J n (E 3 ); 

3. Determine the increment AE 3 from the following 
linear programming problem: 



mm 7, 
AEi 



subject to 



7l J J (E 3 )AE 3 + J (E 3 ) < 7 , 



V T EJ J N ^{E 3 )AE 3 + J N ^(E 3 ) < 7 , 
-0.1 - E 3 <AE 3 <0.1-E 3 . 



4- Let E 3+1 = E 3 + eAE 3 , where e is a small positive 
number controlling the step size; 

5. Repeat Steps (2)-(4) until a desired convergence is 
reached. 

Note that in each iteration we only need to solve a 
linear programming problem, which can be readily cal- 
culated by numerical packages. 



III. OPTIMAL DESIGN RESULTS 

Starting from E(t) = 0.01, after around 12,000 iter- 
ations and 350+ hours computation on a desktop com- 
puter with Intel i5 CPU, we have obtained a satisfactory 
optimal laser pulse as shown in Fig. [3^a). Fig. [3jb), (c), 
(d) show the optimal laser pulse induced NWP evolution 
for \ip g (R,t)\ 2 , \^ u (R,t)\ 2 , and \if> g (R, t)\ 2 + \i> u (R, t)\ 2 , 
respectively. 

Fig. [3]Ja) reveals an interesting physical story. First of 
all, the optimal laser field is a pulse train. After the in- 
ception of H^~ , nearly no electric field is introduced until 
t — 11 fs. In this period, the NWP propagates freely to 
the outer turning point and then turns back, as shown in 




4 3 12 18 20 24 28 32 
Time (fs) 



FIG. 3: (Color online) (a) The designed optimal laser pulse 
train. The evolution of \ip g (R, t)\ 2 (b), \ip u (R,t)\ 2 (c), and 
\tp g (R,t)\ 2 + \ip u (R,t)\ 2 (d) (all in logarithmic scale). 

Fig. G2Jb). The main pulse appears at around t = 11 fs, 
at which time the NWP is moving inward instead of out- 
ward. This is important because if the laser field starts 
interacting with the NWP when it is moving outward, 
part of the wave packets will directly dissociate and the 
subsequent laser pulse has little chance to pull them back 
to the bound states (irij . 

To transfer ip g (R,t) to tpg =0 (R), ip u (R,t) must be me- 
diated. From a closer look at Fig.[3]Jb) and (c), one may 
find that within each oscillation of the electric field, part 
of tp g and ip u are exchanged. The wave function ipu(R, t) 
mainly distributes close to the range R = 3. The quan- 
tity \tp g (R,t)\ 2 + |-0„(i?,i)| 2 gives a smooth evolution of 
the NWP, as depicted in Fig.^d). At the terminal time, 
the NWP has been transferred to the ground vibrational 
state. 

To gain a deeper understanding of the transfer princi- 
ple, we trace the time-dependent probability evolution of 
each individual vibrational state, which can be written 
as 

P»(t)= \(4>g(Wg(R,t))\ 2 , fari/ = 0,-.-,18. (10) 

Fig. [4] shows P v (i) for the first seven vibrational states. 
It is clear that the population of v = increases to 91% 
at the end of the evolution, and the staircase jumps take 
place at the times when the laser pulse is introduced. The 
quick increasing of the ground state population and the 
precipitous dropping of the excited vibrational states in- 
dicate that the laser induced coupling is roughly an anti- 
Stokcs transition: with higher nuclear vibrational 



FIG. 4: (Color online) The time-dependent probabilities for 
the lowest 7 vibrational states. 
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Frequency (a.u.) 

FIG. 5: The frequency spectrum for the designed optimal 
laser pulse train. 



states is excited from lsa g to 2pa u , and then de-excited 
to ls<jg with lower nuclear vibrational states. The proba- 
bility evolution details show more physical scenarios. At 
the beginning v = 1 and v = 2 states have similar prob- 
abilities. After the first laser pulse, the probability of 
v = 1 is halved, whereas the probability of v = 2 does 
not change much. Surprisingly, after the second pulse, 
the probability of v = 1 is doubled and is much larger 
than that of v = 2. After these two pulses, the vibra- 
tional states with v > 2 are already very small, and the 
upcoming third pulse mainly transfers v = 1 to v = 0. 
During the whole process, v = 1 state works as a tempo- 
rary reservoir for storing some population, for ultimately 
maximizing the population of v = 0. 

The frequency spectrum of the obtained laser pulse 
train is shown in Fig. [5] after performing the Fourier 
transform. The main frequency component is around 
0.25. This is consistent with the optimization result in 
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Fig. E](c) , where ip u (R) is mainly excited at the internu- 
clear distance between 2.5 to 3, for the energy gap be- 
tween lstjg and 2pa u at the corresponding intcrnuclear 
distance is around 0.25. The frequency analysis demon- 
strates that the multi-frequency laser pulse train pumps 
and dumps H^~ with different frequency components. 



IV. CONCLUSIONS 

In conclusion, by restricting the laser pulse duration to 
less than 32 fs and confining the electric amplitude within 
0.1, we numerically design an optimal laser pulse train to 
successfully transfer the initial Franck-Condon NWP to 
the ground vibrational state of H^~ with a population 
91%, and the dissociation probability is only 9%. The 
optimal laser pulse train does not act on the NWP until 
the NWP is moving inward. The field-induced Raman 
transition between lsa g and 2pa u transfers the highly 
excited vibrational states to v = directly, or indirectly 
first to v = 1 but finally to v = state. This control 
algorithm can be extended to other molecules. 
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Appendix: Derivation of V B j J„ 

For completeness, we first briefly describe the numer- 
ical procedure to solve the Schrodingcr equation (fTJ). 
We follow the standard split-operator techniques in 
Refs. |3, 0|55|. Let 

E = [Eo, Ei, ■ ■ • , Ek-i], R — [Ro, Ri, ■ ■ ■ , Rn-x\- 

The solution of Eq. (JXJ) can be written as 

i>{R,T f ) =nfjV^ A V(i?,0), (A.l) 



where 



H k = 



T R + V g (R) d gu {R)E k ' 
(R)E k T R + V U (R) 



which is decomposed as 

H k = T + G k , 

where 



T 



T R 
T R 



, G k 



V g {R) d gu (R)E k 
d gu (R)E k V U (R) 



The propagation operator e iHfcAt in Eq. (|A.1[) can be 
calculated by the split-operator method: 



-iH k At 



Substitution of Eq. (|A.2|) into Eq. (|A.1|) yields 

(K-X \ 

i>{R,T f ) = e 



-%T± 



n 

.k=0 



-iG k At-iTAt \ AT 



e zl ^ip(R,0). 

(A.3) 

Here the terms e~ lT ~ and e ~ tTAt can be calcu- 
lated by Fast Fourier Transform (FFT) fBlj ]. Since 
all the four blocks in G k are diagonal matrices, we 

can transform G k into a block diagonal matrix G k 
x) nx 



diag{G°,Gi,-- - ,G^- 1 }, where 



V g (R n ) d gu (R n )E k 
d gu {Rn)E k V u (R n ) 



(A.4) 



Because is symmetric, it can be derived that 

9 k At f V g {R n ) - V u (R n ) 2d gu {R n )E k 
!sm— I — — -a x H ^ — a x , ( , 



(A.5) 



2 V 0k k 

where a x , a z are Pauli matrices, and 

0k = ^(V g (R n ) - V u {R n )) 2 + {2d gu {R n )E k y. (A.6) 

This completes the numerical solution of Eq. (TTJ). 
Now from Eq. ©, we have 



V B J„ 



Rc 



^(Rn) - e ia ^ g (Rn,T f )) 
xV E (e- ia ip g {Rn,T f )) 



(A-7) 



For an element E k in the vector E, it is easy to get 
d 



dE k 



e-^ g (Rn,T f ) 

' ip g (Rn,T f ) - e- ia i> g {Rn,T f ) ° 



(A- 



dE k 



From Eqs. (|A.3[) . we obtain that 



^{R,T f )=e 



d 
dE~k 

x ® c -iG k At c -iTAt 



-iT± 



K-X 

n 

l=k+l 



dE k 



e -iGiAt e -iTAt 



dE k 



k-X 



J e~ iGiAt e _iTA * )e lT ^^(R,0). 



1=0 



6 



From Eqs. (|A.5j) and (|A.6|) . it follows that 

d e -iGfc At = _d_ e -iGZAt£?k_ , _^_ e -^At 



Proceeding further, we obtain 



i At 9 k At fV g ~V u 2d gu E, 

cos ( — 2 — a z H ^ — 

2 2 \ 8k 8k 

. . 8 k At fV g -V u 2d gu E k 
ism — - — „ o z H ttt cr 



2 V 81 



2 ; r> 

k 



and 



dEk 8 k ' 

e Uk = -le 2 ^sin — ^—o x 



dE, 



2 8 k 



Lastly, we need to calculate in Eq. (|A.8[) . Reread- 



ing Eq. ([8]) and defining 



we obtain 



(9a 1 f dq dp 



3E k p 2 + q 2 \ P dE k q dE k 



where 4§- and -M- are none other than the real and 



imaginary parts of the quantity 



• iV-l 



Combining all these equations, we can calculate V £jJ n 
in an explicit manner. 
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